Four-phase patterns in forced oscillatory systems 
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We investigate pattern formation in self-oscillating systems forced by an external periodic perturbation. Ex- 
perimental observations and numerical studies of reaction-diffusion systems and an analysis of an amplitude 
equation are presented. The oscillations in each of these systems entrain to rational multiples of the perturbation 
frequency for certain values of the forcing frequency and amplitude. We focus on the subharmonic resonant 
case where the system locks at one fourth the driving frequency, and four-phase rotating spiral patterns are ob- 
served at low forcing amplitudes. The spiral patterns are studied using an amplitude equation for periodically 
forced oscillating systems. The analysis predicts a bifurcation (with increasing forcing) from rotating four-phase 
spirals to standing two-phase patterns. This bifurcation is also found in periodically forced reaction-diffusion 
equations, the FitzHugh-Nagumo and Brusselator models, even far from the onset of oscillations where the am- 
plitude equation analysis is not strictly valid. In a Belousov-Zhabotinsky chemical system periodically forced 
with light we also observe four-phase rotating spiral wave patterns. However, we have not observed the transi- 
tion to standing two-phase patterns, possibly because with increasing light intensity the reaction kinetics become 
excitable rather than oscillatory. 



I. INTRODUCTION 

Spatially extended systems characterized by the coexis- 
tence of two or more stable states compose a broad class 
of nonequilibrium pattern forming systems. The most com- 
mon multistable systems are those that exhibit bistability, 
(e.g. chemical systems [jj], Q], vertically vibrated granular 
systems [S], and binary fluid convection Spatial pat- 

terns in these systems involve alternating domains of the 
two different stable states, which are separated from each 
other by interfaces or fronts. Bistable systems support a 
variety of patterns from spiral waves to splitting spots and 
labyrinths Q [7], £| g]. In some systems, such as the 
ferrocyanide-iodate-sulfite reaction ||(| and the oxidation 
of carbon monoxide on a platinum surface JTl|], the bistabil- 
ity arises from the nonlinear nature of the system. In other 
systems such as liquid crystals in a rotating magnetic field 
[12 D], [l4[] and periodically forced oscillators §L3jj, the bista- 
bility arises from a broken symmetry. 

Periodically forced oscillatory systems are convenient sys- 
tems for exploring multistability in pattern formation since the 
number of coexisting stable states can be controlled by chang- 
ing the forcing frequency. Applying a periodic force of suffi- 
cient amplitude and at a frequency ujf as ^Wrj> where u)q is the 
oscillation frequency of the unforced system, entrains the sys- 
tem to the forcing frequency. The entrained system has n sta- 
ble states each with the same oscillation frequency but in one 
of n oscillation phases separated by multiples of 2ir/n. We 
refer to the n different phase shifted states as "phase states" of 
the system. 

Recent experiments using the ruthenium-catalyzed 
Belousov-Zhabotinsky reaction forced by periodic illumina- 
tion revealed subharmonic resonance regimes cjf : ujq =2:1, 
3:2, 3:1, 4:1, with two (2:1) , three (3:2,3:1), and four (4:1) 



stable phase states [|[ [161 ], Patterns consisting of alternating 
spatial domains with a phase shift of 7r are observed within 
the 2:1 resonance regime, and three-phase patterns with 
spatial domains phase-shifted by 2ir/3 are observed within 
the 3:1 resonance regime [Q, [l(| |l7|]. The 4:1 resonance is 
more complicated. Adjacent spatial domains may differ in 
phase by either ir or w/2. As a result the asymptotic patterns 
that develop can have four phases, two phases, or a mixture 
of two and four phases. 

In this paper we explore pattern formation in the 4: 1 res- 
onance regimes. In Section II, we describe our experimen- 
tal observations of four phase patterns in the 4: 1 resonance 
band of the forced Belousov-Zhabotinsky reaction. We then 
present an analytical study of the 4:1 resonance [ [T^ , 19 1 in 
Section [ui|. The study is based on a normal form, or ampli- 
tude equation, approach which is strictly valid only close to 
the Hopf bifurcation of the unforced oscillatory system. In or- 
der to test the analytical predictions and to study the behavior 
of forced systems far from the Hopf bifurcation, which is the 
case in the experiments, we conduct numerical studies of two 
reaction-diffusion models (the FitzHugh-Nagumo and Brus- 
selator). We describe the models and results in Section |v|. In 
Section |v| we discuss and compare the analytical and numeri- 
cal results with the experimental observations. 



II. THE PERIODICALLY FORCED 
BELOUSOV-ZHABOTINSKY REACTION 

We use a light-sensitive form of the Belousov-Zhabotinsky 
(BZ) reaction, a chemical reaction system with oscillatory ki- 
netics, to study the 4:1 subharmonic resonance patterns. In 
the experiments, the chemicals of the BZ system diffuse and 
react within a 0.4 mm thick porous membrane. The system is 
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FIG. 1 : The temporal power spectrum of a 4: 1 resonant pattern from 
the BZ experiment for I = 426 W/m 2 . The peak at w/2tt = 0.0154 
Hz is the response at w//4. 



maintained in a non-equilibrium steady state by a continuous 
flow of fresh, well mixed reactant solutions [^0|| on either side 
of the thin membrane where the patterns form. The unforced 
pattern is a rotating spiral wave wave of ruthenium catalyst 
concentration. 

We periodically force the system using spatially homoge- 
neous square wave pulses of light with intensity /, where / 
is the square of the forcing amplitude, and pulse frequency 
u>f (u)f/2n in Hz). We choose the frequency u>f to be ap- 
proximately four times the natural frequency of the unforced 
oscillations. 

To determine the temporal response of a pattern when it is 
periodically perturbed at a particular pair of (I,u)f) parameter 
values we collect a time series of evenly sampled pattern snap- 
shots; a 60 x 60 pixel region of the 640 x 480 pixel image. 
We sample at a rate of approximately 30 frames/oscillation 
and calculate the Fast Fourier Transform for the time series of 
each pixel. The power spectrum of each pixel is determined. 
An average over all pixels provides a power spectrum of a pat- 
tern, as shown in Fig. [l] The 4:1 resonant patterns exhibit a 
dominant peak at uJf/A in the power spectrum. Higher order 
harmonics are also present. 

An example of a 4: 1 resonant pattern observed in the exper- 
iments is shown in Fig. |[ The rotating four-phase spiral wave 
in Fig. ^(a) is the asymptotic state of the system. This image 
is a plot of the phase angle arg(a), where a = a(x, y) is the 
complex Fourier amplitude associated with the oJf/4 mode for 
each pixel (x, y) in the pattern. The four domains (white, light 
gray, dark gray and black) correspond to the four phase states 
with oscillation phases that are shifted by 0, tt/2, it, and 3ir/2 
with respect to the forcing. 

Fig. g(b) is a different representation of the same data. In 
this case the response a at u/ /4 is plotted in the complex plane 
instead of the x — y plane. This representation of the data al- 
lows us to see the distribution of the oscillation amplitude and 
phase at all pixels in the pattern. The four corners of the dia- 
mond shape in Fig. |2|(b) are the four stable phase states. The 
edges of the diamond shape in Fig. ||(b) are formed from pix- 
els at phase-fronts separating adjacent domains. The majority 
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FIG. 2: A rotating four-phase spiral wave observed in the forced BZ 
reaction, (a) A 5.4 x 5.4 mm 2 region of a reactor image showing 
a 4:1 resonant spiral pattern. The white, light gray, dark gray and 
black domains represent the four phase states of the system, (b) A 
plot of the complex Fourier amplitude a at W//4 for each pixel of 
the pattern. The forcing intensity is 7= 426 W/m 2 and the forcing 
frequency is cuf /2n — 0.062 Hz. The data were filtered to isolate 
the response at U//4 from the higher harmonics. 
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FIG. 3: A histogram showing the distribution of phase angles in the 
pattern in Fig. 0. The four peaks indicate the high density of points 
in each of the four phase states. 



of pixels in the pattern are in one of the four corner states as 
the histogram of phase angles in Fig. ^illustrates. 

Traveling four-phase patterns exist over the entire dynamic 
range of forcing intensity I in the 4: 1 resonance region. The 
range of forcing intensity is limited by a /-dependence of the 
reaction kinetics. As / is increased, the reaction kinetics shifts 
from oscillatory to excitable. 



III. AN AMPLITUDE EQUATION FOR FORCED 
OSCILLATORY SYSTEMS 

We study the experimental observations shown in the pre- 
vious section using a normal form equation for the amplitude 
of the u)f/4 mode. Consider first an oscillatory system re- 
sponding to the forcing at lOf/n where n is integer. We as- 
sume the system is near the onset of oscillations, i.e. close to 
a Hopf bifurcation. The set of dynamical fields u describing 
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the spatio-temporal state of the system can be written as 

u = u Aexp (iujft/n) + c.c. + . . . , (1) 

where Uo is constant, A is a slowly varying complex ampli- 
tude, and the ellipses denote other resonances with smaller 
contributions. The slow space and time evolution of the ampli- 
tude A is described by the forced complex Ginzburg-Landau 
equation, 

(jj, + iv)A + (1 + ia)A zz - (1 - i(3)\A\ 2 A 



A T 



InA 



*(n-l) 



(2) 



where [i is the distance from the Hopf bifurcation, v is the 
detuning from the exact resonance, and j n is the forcing am- 
plitude. 

For the special case n — 4 (the 4: 1 resonance) we can elim- 
inate the parameter fj, by rescaling time, space, and amplitude 
as t = /it, x = \J fi/2 z and B — A^/Ji to obtain 

(1 + iis )B + hi + ia)B xx - (1 - i/3)\B\ 2 B 



B, 



+ lB*% (3) 

where i/q = v/fi. Equation ([|) also applies to the 4:3 sub- 
harmonic resonance. This follows from symmetry consider- 
ations: the system is symmetric to discrete time translations 
t^t + ^jj=t+7^. The amplitude equation must then be 



invariant under the transformation B 



Bcxp(3ni/2). The 
is 



only forcing term satisfying this requirement to cubic order 



A. Phase states and phase fronts 

Constant solutions of Eq. ^ indicate that the system is en- 
trained to the forcing. There are four stable constant solutions 
to Eq. (Q), each with the same amplitude but with different 
phases, axg(B), which correspond to the four stable phase 
states. Simple expressions for these solutions and exact forms 
for the front solutions connecting them in space are obtained 
from the gradient version of Eq. where vq = a — f3 — 0: 



Bt 



B + -B T 



\B\ 2 B 



■ 7 B* 



(4) 



The stable phase states (constant solutions) of Eq. (|4|) for < 
7 < 1 are (B 1: B 2 , B 3 , _B 4 ) = (A, iX, - A, -iA) where A = 
— 7- They are represented as solid circles in Fig. 4. 
Front solutions connecting pairs of these states are of two 
types, fronts between states separated in phase by tt and fronts 
between states separated in phase by ir/2 (hereafter 7r-fronts 
and 7r/2-fronts). The 7r-front solutions are 



i?3_>i = -Bitanha;. 
-B4_>2 = -E>2 tanh x . 

For the particular parameter value 7 = 
have the simple forms 



(5) 

1/3 the 7r/2-fronts 



B 2 ^i = ^J|"[l + *+(l-i)tanhar] , 
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FIG. 4: The four phase states (black dots) connected by phase-fronts 
in the forced complex Ginzburg-Landau equation (Q). Two types of 
fronts between phase states are possible; the solid lines are 7r-fronts 
and the dashed lines are 7r/2-fronts. 
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Additional front solutions follow from the invariance of 
Eq. (Q) under reflection, x — > —x. 

Figure 4 shows these front solutions (parametrized by the 
spatial coordinate x) in the complex B plane. For example, 
the 7r-front i?3^i is represented by the solid line connecting 
the state £> 3 to the state B L as x increases from — 00 to +00. 
The 7r/2-front B 2 ^i is represented by the dashed line con- 
necting the state B 2 to the state B\. 

In the special case of the gradient system (Q) all front solu- 
tions are stationary. The more general case with nongradient 
terms in Eq. (^]) can be studied by perturbation theory when 
1/0, a and f3 are small []19|], The results of this analysis show 
that the 7r/2-fronts become propagating fronts while the ir- 
fronts remain stationary. 

Figure || shows a rotating four-phase spiral wave from a 
numerical solution of the two-dimensional version [ pT| ] of 
Eq. (^|). The phase diagram in the complex B plane, shown 
in Fig. ||(b), has four 7r/2-fronts: -Bi_>4, £?4^3, .63^2 and 
B 2 ^\. The amplitude B corresponds to the complex Fourier 
amplitude a measured in the experiment; the four-phase spiral 
pattern in Fig. [|and the corresponding diamond-shape in the 
complex plane are predicted by the amplitude equation. 



B. A phase-front instability 

The existence of the stationary 7r-front solutions suggests 
that standing two-phase patterns similar to those found un- 
der 2:1 resonant conditions [^, |l6[] may be observed in the 
4:1 resonant case provided the 7r-fronts are stable. Standing 
two-phase patterns have not been observed in experiments in 
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FIG. 5: A rotating four-phase spiral wave in the forced complex 
Ginzburg-Landau equation, (a) arg(_B(x, y)) in the x — y plane, (b) 
B(x,y) in the complex plane. Parameters: 7 = 0.6, v$ = 0.1, 

a = (3 = 0. 



the 4: 1 resonance band so the stability of 7r-fronts becomes a 
question. Stability conditions for 7r-front solutions were stud- 
ied in Refs. [[18, 19 1. The results are described below. 

Consider the pair of 7r/2-fronts shown in Fig. ||(a). They 
are separated by a distance 2\ and connect the phase states 
£?3 and B\. For 7 1/3, the solutions (Q) are good ap- 
proximations to 7r/2-front solutions. The pair of fronts can 
be represented as 



B(x;Cx) « B 3 ^ 2 (x-(+x)+B 2 ^ 1 (x-(- X )-i\, (7) 

where C, is their mean position. For large separation distances 
(X >> 1) B w i?3^2 when x « C — X ar, d -S ~ -B2_>i 
when x w C + X an d Eq. (0) represents a pair of isolated 7r/2- 
fronts. When the distance between the pair decreases to zero 
(X — > 0), then B w Ba^i and Eq. (^|) approaches a 7r-front 
solution. 

The stability of 7r-fronts is determined by the interaction 
between a pair of 7r/2-fronts. Stable 7r-fronts are the result 
of an attractive 7r/2-front interaction; the 7r/2-fronts attract 
each other and the distance between them decreases to zero. A 
repulsive interaction implies unstable 7r-fronts. The potential 
V(x) that governs this interaction, 



dV 
dx ' 



(8) 



is shown in Fig. ^(b) for various 7 values. The potential has 
a single maximum for 7 < -f c = 1/3 which represents a re- 
pulsive interaction between 7r/2-fronts and the instability of 
7r-fronts. It has a single minimum for 7 > j c which indicates 
the attractive interaction between 7r/2-fronts and the resulting 
stability of 7r-fronts. At j c the potential is flat, V — 0, for 
all x values. At this parameter value, pairs of 7r/2-fronts do 
not interact and there is a continuous family of front pair solu- 
tions with arbitrary separation distances, 2%, in Eq. ([7]). This 
degeneracy of solutions at the critical point 7 = j c is removed 
by adding higher order t erms to the amplitude equation, as we 



discuss in Section [II D 
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To summarize, stationary 7r-front solutions of Eq. (|3|) are 
stable for forcing amplitudes 7 > -f c = 1/3. When 7 is 
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FIG. 6: (a) The phase, arg(B) of a pair of 7r/2-fronts, i?3_,2, and 
The distance between the two fronts is defined to be 2\- 
(b) The potential V(x) describing the interaction between two 7r/2- 
fronts. For 7 > y c the extremum at x — is a minimum and \ 
converges to 0. For 7 < j c the extremum is a maximum and \ 
diverges to ±00. At 7 = 7c the potential is flat and there is no 
interaction between 7r/2-fronts. 



decreased past j c , 7r-fronts lose stability and split into pairs 
of propagating 7r/2-fronts. The splitting process is shown in 
Fig. [7] where the B^\ 7r-front evolves into the pair of stable 
traveling 7r/2-fronts, B^ 2 and B 2 ^\ when 7 < j c . The par- 
ity symmetry x —> — X makes evolution toward the pair Si_>4 
and £?4^3 equally likely. The splitting occurs for forcing am- 
plitudes arbitrarily close to j c , although in that case the time 
scale of this process becomes very long. 



C. Effects of the phase-front instability on pattern formation 

The stability of stationary 7r-fronts for 7 > j c suggests the 
predominance of standing two-phase patterns. These patterns 
involve alternating domains with oscillation phases shifted 
by 7r with respect to one another. Domains shifted by ir/2 
may exist as transients; the interactions between 7r-fronts and 
7r/2-fronts always produce 7r/2-fronts which are stable but 
attract one another and coincide to form stationary 7r-fronts. 
Since the 7r/2-fronts are traveling these transients are rela- 
tively short. For 7 < j c the interactions between the tt/2- 
fronts are repulsive. The 7r-fronts are unstable and split into 
pairs of traveling 7r/2-fronts. As a result, traveling waves with 
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FIG. 7: An example of the phase-front instability in one space di- 
mension. Left: The space-time plot (solutions of Eq. (Q)) shows the 
splitting of an unstable 7r-front into a pair of traveling 7r/2-fronts. 
The 7r/2-front pairs enclose the dark gray domain that has an os- 
cillation phase shifted by n/2 with respect to the black and light 
gray domains. Right: Snapshots at times t = 0, t = 100, and 
t = 300, showing the instability in the complex B plane. Param- 
eters in Eq. (|): u = 0.02, 7 = 0.3, a = (3 = 0. 



all four phase-states are the asymptotic pattern. 

A typical two-dimensional traveling pattern involving all 
four phases is the four-phase spiral wave shown in Fig. |] or 
in Fig. U Figure shows the effect of the phase-front in- 
stability on a four-phase spiral wave. The initial spiral wave 
(Fig. |8|(a)) was obtained by solving a two-dimensional version 
of Eq. ([3]) for 7 < j c . The following three frames (Fig. |^(b)- 
(d)) are snapshots showing the evolution of the initial four- 
phase spiral wave into a standing two-phase pattern after 7 is 
increased above j c . The evolution begins at the spiral core 
where the attractive interactions between pairs of 7r/2-fronts 
are the strongest. The coalescence of 7r/2-fronts leaves behind 
a stationary 7r-front which grows in length until no 7r/2-fronts 
are left, as is evident by the single line in the complex B plane 
shown in Fig. ^|(d). 



D. Higher order terms in the amplitude equation 

From the analysis of Eq. (^) we have shown that two-phase 
patterns must be standing and four-phase patterns must be 
traveling. The analysis of the equation with higher order con- 
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FIG. 8: Numerical solution of a two-dimensional version of Eq. (|3j) 
showing the evolution of a rotating four-phase spiral-wave into a 
standing two-phase pattern when 7 is increased above y c . The frames 
on the left show arg(B) in the x — y plane. The frames on the right 
show the complex B plane, (a) The initial four-phase spiral wave 
(computed with 7 < y c ). (b) The spiral core, a 4-point vertex, splits 
into two 3-point vertices connected by a 7r-front. (c) A two-phase 
pattern develops as the 3-point vertices further separate, (d) The fi- 
nal standing two-phase pattern. Parameters: 7 = 0.6, Vo = 0.1, 
a = f3 = 0, 7c w 1/3. 



tributions suggests the possible existence of a small 7 range, 
of order /1 <C 1, surrounding j c where slowly traveling two- 
phase patterns exist. 

The higher order contributions to Eq. (Q), such as |-B| 4 5, 
or \B\ 2 B XX , lift the de generacy of the instability. Figure H 
shows two possible scenarios for the front interaction poten- 
tial V when higher order contributions to Eq. ([|) are included 
(both scenarios lift the degeneracy of the phase-front instabil- 
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ity). In one case, shown in Fig. |](a), the stationary 7r-front 
loses stability to a pair of counter-propagating 7r-fronts in a 
pitchfork bifurcation which leads to double-minimum poten- 
tial. This scenario is a nonequilibrium Ising-Bloch pitchfork 
bifurcation of 7r-fronts like the one found in the 2: 1 resonance 
case §|] and in other bistable systems ||, |4[ |||. It 
leads to slow traveling two-phase patterns in the range where 
7 is near j c . In the scenario shown in Fig. |](b), the stationary 
7r-front loses stability via a subcritical bifurcation which leads 
to double-maximum potential. In this case there is a range of 
stable 7r-fronts coexisting with pairs of separated 7r/2-fronts. 
This allows the possibility of patterns containing both 7r-fronts 
and 7r/2-fronts. Beyond this range the potential has a single 
maximum and 7r-fronts split into pairs of 7r/2-fronts. Both 
scenarios persist over a range of 7 of order /1, the distance 
from the Hopf bifurcation. 
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FIG. 9: The degeneracy of the potential V (x) at 7 = 7 C is broken by 
adding higher order terms to (^). In the intermediate range of 7 « j c 
two scenarios are possible as 7 is decreased through the bifurcation: 
(a) The \ — solution loses stability in a pitchfork bifurcation at 
7 C to a pair of solutions that move to ±00. (b) The \ — solution 
remains stable while the \ — i 00 solutions acquire stability and 
lose stability only below j c . In both cases the deformations from a 
single minimum to a single maximum occur within a small range of 
7 of order p< 1. 



IV. NUMERICAL SOLUTIONS OF PERIODICALLY 
FORCED REACTION-DIFFUSION MODELS 

The amplitude equation analysis predicts the existence of a 
phase-front instability near the Hopf bifurcation and hints at 
possible modifications of the instability as the distance from 
the Hopf bifurcation is increased. Our objectives in this sec- 



tion are to test the existence of the instability in reaction- 
diffusion models and to use the models to examine how the 
instability is modified far from the Hopf bifurcation. 



A. The FitzHugh-Nagumo model 

We study a periodically forced version of the FitzHugh- 
Nagumo equations 



u t = u — (1 + Tcosuift)i 
vt = e(u — aiv) + 5V 2 v . 



v + W 2 u. 



(9) 



The unforced model is obtained by setting V = 0. The 
uniform state (u, v) = (0, 0) undergoes a Hopf bifurcation 
as e is decreased past e c — l/oi. The Hopf frequency is 
= V £ c — 1 and the distance from the Hopf bifurcation is 
measured by \i = (e c — e)/e c . 

We compute the numerical solutions of Eq. (^) in the 4: 1 
resonance band (cof ps 4cj#) and close to the Hopf bifurca- 
tion (fx <C 1). Close to the Hopf bifurcation the amplitude 
equation analysis applies. We expect to find a critical value of 
the forcing amplitude r c corresponding to the phase-front in- 
stability point 7 C in the amplitude equation. For the FitzHugh- 
Nagumo equations this r c will, in general, depend on the pa- 
rameters e, S, a%, and uoj. In the following we fix a-y = 1/2, 
5 = 0, u>t = 4 and only vary e (the parameter that controls the 
distance \l to the Hopf bifurcation) and the forcing amplitude 

r. 

Close to the Hopf bifurcation we find stable stationary 7r- 
fronts for forcing amplitudes T > T c . Below T c , stationary ir- 
fronts are unstable and split into pairs of 7r/2-fronts. Figure |lC| 
illustrates this in numerical solution of a one-dimensional ver- 
sion of Eq. (^). An stable 7r-front pattern is generated from 
random initial conditions with T > T c . At t = T is de- 
creased below T c ; the 7r-front becomes unstable and splits 
into a pair of traveling 7r/2-fronts. 

The numerically computed r c for the solution in Fig. nQis 
r c ps 2.15. Since T c is a function of the parameters in Eq~^|), 
we define a new parameter r\ — (r c — T)/r c that measures 
the distance from the phase-front instability point. In Fig. [To], 
77 w 0.012 indicating that we are just beyond the critical point. 

Farther from the Hopf bifurcation we find that the phase- 
front instability still exists. Figure [IT] shows the the evolution 
of an initial unstable stationary 7r-front with parameters cho- 
sen so the system is far from the Hopf bifurcation but at the 
same distance, r\ « 0.012, from the phase-front instability. 
The asymptotic solution is a slowly propagating 7r-front, in 
contrast to a pair of separated 7r/2-fronts that develop close 
to the Hopf bifurcation (see Fig. [R)|). The range of forcing 
amplitudes near T c over which these traveling 7r-fronts exist 
increases with /1. At smaller forcing amplitudes, below the 
range of traveling 7r-fronts, 7r-fronts split into pairs of tt/2- 
fronts and four phase traveling patterns prevail. 

In two dimensions the typical traveling wave pattern for 
r < T c is a rotating four-phase spiral wave. Figure |l2|(a) 
shows a stable four-phase spiral wave generated from random 
initial conditions. Using this spiral as an initial condition, we 
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FIG. 10: The phase-front instability in the 4:1 resonance of the 
forced FitzHugh-Nagumo model close to the Hopf bifurcation. Left: 
a space-time plot of arg(a) where a is the complex Fourier coeffi- 
cient of the 4:1 response (equivalent to A in Eq. (Q)). At t = the 
forcing amplitude was decreased below T c . The initial standing n- 
front becomes unstable and splits into a pair of traveling 7r/2-fronts. 
The 7r/2-fronts separate the black, dark gray, and light gray domains 
where the oscillation phase is shifted successively by tt/2. Right: 
The same data depicted in the complex a plane at three successive 
times, t = 0, t = 5601/ , and t = 4160T/ where Tf = 2-K/u f . 
(a) The initial standing 7r-front is unstable, (b) The front develops 
an intermediate phase, (c) Two tt/2 fronts are formed. Parameters: 
ai = 0.5, e = 1.95, 8 = 0, T = 2.0, u f = 4.0, and fj, = 0.025. 
The phase-front instability point is T c ~ 2.15 and r\ ~ 0.012. 
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FIG. 11: The phase-front instability in the 4:1 resonance of the 
forced FitzHugh-Nagumo model far from the Hopf bifurcation. Left: 
a space-time plot of arg(a). At t = the forcing amplitude was de- 
creased below r c . The initial standing 7r-front is unstable and and 
starts traveling to the right. In this case, no intermediate phase de- 
velops. Right: The same data depicted in the complex a plane at 
three successive times, t = 0, t = 30001/, and t = 60001/ where 
Tf = 2ir/ujf. (a) The initial standing 7r-front. (b) The standing 
7r-front is unstable and begins to travel, (c) The asymptotic pattern 
is a traveling 7r-front. Parameters: a 1 — 0.5, e — 1.5, 5 = 0, 
r = 1.585, u>f = 4.0, and [i — 0.25. The phase-front instability 
point is T c w 1.605 and r\ » 0.012. 



B. The Brusselator model 



increase T above T c and the system evolves into a two-phase 
standing pattern. Figures [l^(b)-(d) show the transition. Since 
the 7r/2 fronts are attracting the spiral is unstable and two of 
the four phase domains shrink until a standing two-phase pat- 
tern remains. 

The numerical solutions of the forced FitzHugh-Nagumo 
equations support the predictions of the amplitude equation 
analysis. Close to the Hopf bifurcation, the phase-front in- 
stability is found (compare Fig. ^ with Fig. [R^ and Fig. |8] with 
Fig. [j~2|). Far from the Hopf bifurcation the instability persists. 
The effects of higher order terms in the amplitude equation 
are valid even far from the Hopf bifurcation (jjl = 0.25); the 
phase-front instability near the Hopf bifurcation (as // — > 0) 
turns into an Ising-Bloch pitchfork bifurcation. Stationary ir- 
fronts bifurcate to traveling 7r-fronts and not 7r/2-fronts. 



We tested the transition from four-phase traveling waves 
to two-phase standing waves using another reaction-diffusion 
model, the forced Brusselator, 

u t = c-du+ [1 + T cos uj f t]u 2 v + V 2 u, (10) 
v t = du — u 2 v + (5V 2 u . 

The unforced Brusselator, obtained by setting T = 0, has a 
stationary uniform state (u, v) = (c, d/c) which undergoes 
a Hopf bifurcation as d is increased past d c = 1 + c 2 . The 
Hopf frequency is u>h = c and the distance from the Hopf 
bifurcation is measured by fj, = (d — d c )/d c . 

We studied Eq. (|o|) in the 4: 1 resonance band using a nu- 
merical partial differential equation solver [^(| We found 
that below a critical forcing amplitude T c the solutions are ro- 
tating four-phase spiral waves consisting of 7r/2-fronts (see 
Fig. [T^(a)). The four-phase spiral wave was generated by one 
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FIG. 12: Numerical solution of the forced FitzHugh-Nagumo equa- 
tions (^) in 4:1 resonance shown at four successive times t = 0, 
t = 11600T/,t = 13600T/, andt = 15600T/ where Tf = 2ir/cj f . 
The frames on the left show arg(a) in the x — y plane. The frames 
on the right show the complex a plane, (a) The initial spiral wave 
of four phases separated by 7r/2-fronts is computed with Y < F c . 
(b) When V is increased above r c two pairs of 7r/2-fronts begin to 
attract one another, (c) As the 7r/2-fronts attract they collapse into 
a stationary 7r-front which grows in length, (d) The final pattern is 
two phase domains separated by a stationary 7r-front. Parameters: 
ai = 0.5, e = 1.5, 5 = 0, T = 2.5, and uj f = 4.0. 

of two following initial conditions: a spiral wave computed 
from the unforced (T = 0) Brusselator equations, or the linear 
functions 

u{x,y)= y/L, 0<y<L, 
v(x,y)= -2x/L + A < x < L , 



Above r c pairs of it/ 2 fronts attract each other and the core 
of the spiral evolves into an expanding 7r-front. Figures [l3|(b)- 
(d) illustrate this process. When the ir/2 fronts disappear, the 
resulting asymptotic pattern is two states separated by a sta- 
tionary 7r-front. The transition from a four-phase spiral wave 
to a two-phase stationary pattern, as in the amplitude equation 
model and the FitzHugh-Nagumo model, indicates the exis- 
tence of the phase-front instability in the Brusslator model. 
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FIG. 13: Numerical solutions of the forced Brusselator model ( |10| ) 
showing snapshots at t = 0, t = 748T), t = 10007/ , and t = 
55447 1 / where Tf = 2n/uj>f. The rotating four-phase spiral wave, 
computed with F < F c (T — 0.11), transforms into a standing two- 
phase pattern after T is increased past r c (T = 0.13). The frames 
in the left column show arg(a) in the x — y plane where a is the 
complex Fourier coefficient of the 4: 1 mode. The right column shows 
the same data in the complex a plane. Parameters: c = 0.5, d — 1.5, 
5 = 5.0, ujf = 1.69, and fi — 0.20. The numerical solution grid 
was 128 x 128 points. 



where L — 632.5. 
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V. CONCLUSIONS 

We studied 4: 1 resonant patterns in Belousov-Zhabotinsky 
chemical experiments, in an amplitude equation for forced 
oscillatory systems (the forced complex Ginzburg-Landau 
equation), and in forced FitzHugh-Nagumo and Brusselator 
reaction-diffusion models. At low forcing amplitudes all of 
these systems exhibit traveling four-phase patterns. 

An analysis of a forced complex Ginzburg-Landau equa- 
tion, derivable from periodically forced reaction-diffusion 
systems near a Hopf bifurcation, predicts traveling four-phase 
patterns at low forcing amplitude and standing two-phase pat- 
terns at high forcing amplitude. The transition mechanism 
between these two patterns is a degenerate phase-front insta- 
bility where a stationary 7r-front splits into a pair of traveling 
7r/2-fronts. We derived an interaction potential between w/2- 
fronts that describes the instability as a change from repulsive 
to attractive 7r/2-front interactions. We investigated the be- 
havior of the instability near the critical point where higher 
order terms in the amplitude equation become important. We 
found that these terms lift the degeneracy of the instability 
and introduce a narrow intermediate regime. In this regime 
we found both slowly traveling 7r-fronts and the coexistence 
of stable stationary 7r-fronts and repelling pairs of 7r/2-fronts. 

We further investigated this phase-front instability using 
the FitzHugh-Nagumo and the Brusselator reaction-diffusion 
models. These models exhibit the instability even far from the 
Hopf bifurcation where the amplitude equation is not known 
to be valid. Near the Hopf bifurcation the instability, at T c , 
separates patterns of stationary 7r-fronts from patterns of trav- 
eling 7r / 2-fronts. In two dimensions, a rotating four-phase spi- 
ral wave evolves into a two-phase standing pattern when T is 



increased past T c . In the FitzHugh-Nagumo model we found, 
far from the Hopf bifurcation, an intermediate range near T c 
where traveling it front patterns were observed. These numer- 
ical results are in full agreement with the theoretical predic- 
tions based on the amplitude equation. 

The standing two-phase patterns found in the amplitude 
equation and in the FitzHugh-Nagumo and Brusselator mod- 
els were not observed in the experiments, which were con- 
ducted far from the Hopf bifurcation. However, the existence 
of the phase-front instability far from the Hopf bifurcation was 
found in the numerical studies of the FitzHugh-Nagumo and 
Brusselator models. We conclude that the large distance from 
the Hopf bifurcation does not explain the absence of standing 
two-phase patterns in the experiments. A more likely expla- 
nation is the limited dynamic range of the forcing amplitude 
in the experiments. Experiments show that the dynamics of 
the BZ reaction are 7-dependent; as the forcing amplitude is 
increased, the dynamics undergo a transition from oscillatory 
to excitable kinetics. The excitable kinetics are not described 
by the amplitude equation or by the reaction-diffusion models 
in the parameter ranges we studied. 
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